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ABSTRACT 



£N| ■ Aims. We present a weak lensing search of galaxy clusters in the 4 deg 2 of the CFHT Legacy Survey Deep. This work aims at building a 
mass-selected sample of clusters with well controlled selection effects. This present survey is a preliminary step toward a full implementation 
in the forthcoming 170 deg 2 of the CFHTLS Wide survey. 

Methods. We use the deep i' band images observed under subarcsecond seeing conditions to perform weak lensing mass reconstructions 

in ■ 

and to identify high convergence peaks. Thanks to the availability of deep u'g'r'i'z' exposures, sources are selected from their photometric 
redshifts in the weak lensing analysis. We also use lensing tomography to derive an estimate of the lens redshift. After considering the raw 
statistics of peaks we check whether they can be associated to a clear optical counterpart or to published X-ray selected clusters. 
Results. Among the 14 peaks found above a signal-to-noise detection threshold v = 3.5, nine are secure detections with estimated redshift 
0. 1 < zi < 0.7 and a velocity dispersion 450 < <x,, < 700kms~'. This low mass range is accessible thanks to the high density of background 
sources. Considering the intersection between the shear- selected clusters and XMM-LSS X-ray clusters in the Dl field, we observe that the 
ICM gas in these low-mass clusters (T x ~ 1 — 2 keV) is not hotter than the temperature inferred from shear, this trend being different for 
published massive clusters. A more extended weak lensing survey, with higher statistics of mass structures will be a promising way to bypass 
several of the problems related to standard detection methods based on the complex physics of baryons. 
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■ 1. Introduction have an indirect access to cluster masses and rather measure 
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In the context of hierarchical structure formation within the 



quantities like the SZ decrement, X-ray or optical/NIR lumi- 
nosities, the line-of-sight velocity dispersion of the cluster- 
Lambda Cold Dark Matter paradigm (ACDM), clusters of , , . . , ... , , . . . 

r b v ' member galaxies or the temperature of the hot gas in the mtra- 
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2001; Borgani & Guzzo 2001). 



quantities like mass and because the process of detecting clus- 
ters of galaxies with such indirect methods might suffer var- 

Therefore the main physical parameter for a cluster of ioufj kinds of selection effectfj Therefore any attempt t0 use 

galaxies is its total mass. However most of observations only dusters Qf galaxiefj ag effident cosmological probes cannot af . 
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Since the early 90's several groups have reported the detec- 
tion of a weak lensing signal around massive clusters of galax- 
ies. However the broad range of observational configurations 
(field of view, depth, seeing, ground- or space-based images, 
etc) makes difficult a direct comparison of published results. 
For reviews see Mellier (1999) and Bartelmann & Schneider 
(2001). Progresses have been made in this direction with weak 
lensing studies of sizable samples of optically or X-ray selected 
clusters (Dahle et al. 2002; Cypriano et al. 2004; Clowe et al. 
2006; Bardeau et al. 2006). Different mass estimates globally 
agree although outliers perturb a simple relation between X-ray 
(or dynamical) and lensing mass estimates (e.g. Allen 1998; 
Wu 2000; Arabadjis et al. 2004). This suggests that dynamical 
activity is still important for massive halos and that aspheric- 
ity and projection effects may complicate both weak lensing 
and other mass estimates (Metzler et al. 2001; Hoekstra 2003; 
Clowe et al. 2004; de Putter & White 2005; Gavazzi 2005; De 
Filippis et al. 2005). 

In parallel, the idea of direct detection of galaxy clusters by 
their weak lensing signal starts to emerge. By measuring the co- 
herent stretching of distant galaxies by intervening structures, 
one is able to infer the projected density field (i.e. the so-called 
convergence). Hence high convergence peaks may be identified 
as massive clusters of galaxies. This is the idea of a direct weak 
lensing cluster survey (hereafter WLCS), aimed at building a 
mass-selected cluster sample directly comparable to CDM the- 
ory (through N-body cosmological simulations). On the theo- 
retical side, pioneering analytical predictions based on the mass 
function of halos have been proposed (Schneider 1996; Kruse 
& Schneider 1999), but they were not able to properly account 
for projection effects (Reblinsky & Bartelmann 1999). In ad- 
dition Bartelmann et al. (2001) showed that WLCSs are very 
sensitive to the details of the clusters density profile. More re- 
cently, ray-tracing into N-body cosmological simulations have 
been used to properly address the critical issue of projections 
and clusters' asphericity (White et al. 2002; Padmanabhan et al. 
2003; Hamana et al. 2004; Tang & Fan 2005) and the way 
to reduce the effect of noise on cluster detections through an 
optimised data filtering procedure (Hennawi & Spergel 2005; 
Maturi et al. 2005; Starck et al. 2006). A step forward will 
also probably be made with simplified analytical models of 
the convergence one-point PDF (Taruya et al. 2002; Das & 
Ostriker 2006). Under standard observational conditions, these 
works predict that at z ~ 0.2 clusters more massive than 
M ~ 5 x 10 13 M G can be recovered with a signal-to-noise ra- 
tio v = 3. This limit drops to M ~ 2 x 1O 14 M at z ~ 0.7 . 
Therefore the main targets of WLCSs are massive clusters of 
galaxies. 

From the point of view of observations, we can mention 
a few serendipitous detections of galaxies clusters via weak 
gravitational lensing (e.g. Schirmer et al. 2003, 2004; Dahle 
et al. 2003). A few examples have also been found to show 
up through weak lensing techniques without any clear optical 
counterpart and gave support for the existence for the so-called 
"dark clumps" (Umetsu & Futamase 2000; Miralles et al. 
2002). The practical implementation of a systematic WLCS is 
however very new. Miyazaki et al. (2002) studied an area of 
2.1 deg 2 with Suprime-Cam on Subaru telescope under excel- 



lent seeing conditions. They report an excess of 4.9 ± 2.3 con- 
vergence peaks with signal-to-noise v > 5. Hetterscheidt et al. 
(2005) report the detection of ~ 5 cluster candidates over a set 
of 50 disconnected VLT/FORS deep images covering an effec- 
tive area of 0.64 deg 2 while Wittman et al. (2006) present pre- 
liminary results for the first 8.6 deg 2 of the Deep Lens Survey 
(eight detections). Haiman et al. (2004) also make interesting 
predictions for future WLCS applications in the LSST survey. 

In this paper we present a weak lensing analysis of the Deep 
CFHT Legacy Survey 1 covering 4x1 deg 2 in five optical bands 
(u*g'r'i'z') under subarcsec seeing condition as a pilot analysis 
for the ongoing Wide Survey which will cover 170 deg 2 . The 
present work proposes to carry out weak lensing mass recon- 
structions in the Deep fields and focus on high convergence 
peaks in order to shed light on WLCS capabilities. The rela- 
tively high sample variance of the Deep images prevents any 
cosmological application of WLCSs but the great depth and 
amount of photometry make us with an excellent laboratory 
for a future application to the Wide data. 

The paper is organised as follows. In Sect. 2 we briefly re- 
view the basics of weak gravitational lensing. In Sect. 3 we 
present the data at hand, the specific treatment required for 
weak lensing signal extraction and photometric redshifts. We 
show mass reconstructions (i.e. convergence maps) in Sect. 

4 inferred from the coherent shear field imprinted on distant 
background sources. We also measure the statistics of high con- 
vergence peaks whereas we focus on their properties in Sect. 

5 by studying the associated optical and X-ray counterparts 
(when available). We discuss our main achievements and con- 
clude in Sect. 6. 

In the following we assume the "concordance model " cos- 
mological background with Ho = 70 km s Mpc~', O m = 0.3 
and O a = 0.7 . All magnitudes are expressed in the AB photo- 
metric system. 

2. Basic lensing equations 

In this section we briefly summarise the necessary background 
of gravitational lensing and especially the weak lensing regime 
which concerns the present analysis. We refer the reader to the 
reviews of Mellier (1999) and Bartelmann & Schneider (2001) 
for more detailed accounts. 

The fundamental quantity for gravitational lensing is the 
lens Newtonian potential ifr(0) at angular position which is 
related to the surface mass density 1.(6) projected onto the lens 
plane through 

4GDoiAs f d 2 (ri:ie > )]n]0 _ > l! (1) 

c z D os J 

where D i, D os and Dj s are angular distances to the lens, to the 
source and between the lens and the source respectively. The 
deflection angle a = Viff relates a point in the source plane /? 
to its image(s) in the image plane through the lens equation 
P = - a(0). The local relation between j8 and is the Jacobian 
matrix a^ = 8(3,1 89 j 

1 http : //www . cfht . hawaii . edu/Science/CFHTLS/ 
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The convergence k(6) = E(0)/E crit is directly related to the sur- 
face mass density via the critical density 

c 2 As 

Scrit = ^ " . (3) 



4ttG D iAs 
and satisfies the Poisson equation 

Aiff = \jj n + iff 22 = 2k 
_ l 



(4) 

The 2-component shear is y = j\ + iy 2 = \(^,n - <A,22) + ^,12 
in complex notation. An elliptical object in the image plane is 
characterised by its complex ellipticity e defined from the sec- 
ond moments tensor Q t j of the surface brightness distribution 

as 

Qu - Q22 + 2iQ n 



A/2 



(5) 



211 + 222+2 (QnQ22-Q 2 u ) 

This observed ellipticity can be related to the intrinsic elliptic- 
ity of the source e s by: 

e s + g 



1 + g*e s 
e* s +g* 



,for \g\<l 



,for \g\ > 1 



(6) 



where g = y/(l - k) is the so-called reduced shear and * de- 
notes complex conjugation (Seitz & Schneider 1997; Geiger & 
Schneider 1998). In the weak lensing regime (g ~ y «c 1), Eq. 
(6) reduces to e — e s + y. Provided the random orientation of 
sources reduces the averaged source ellipticity to zero, e pro- 
vides an unbiased estimate for the shear y. The noise associated 
to this estimator is due to the scatter in the intrinsic ellipticity 
of sources with a typical value cr e ~ 0.3 per component. 

In the equations above we can isolate a geometric term 
which linearly scales the lensing quantities k, iff, and y and 
only depends on the distance ratio As /Ads- We thus can write 
y = w(zi> Zs)Too (and so forth for k and iff) with w(z\,z s ) = 
As/As®(z s - Zi) and 0(x) is the Heavyside step function. If 
sources are not confined in a thin plane, we account for the dis- 
tribution in redshift by defining an ensemble average distance 
factor W(zi) such that: 



W( Z1 ) = (w(zuzs)) zs = 



f 

•Jz\ 



dz s nUs) — 

^os 



f 

JO 



(7) 



dz s n(z s ) 



If we now consider a broad distribution of mass interven- 
ing between sources and the observer, we can use the Born ap- 
proximation for the propagation of light in a clumpy Universe 
to infer the effective convergence experienced by light bundles: 

^p(l,0)-p(l) 



*eff(0) 



'crit(0 
2 rXH 



dl 



dxW(x)fk(x) 



s(fk(x)e,x) 

a(x) 



(8) 



with I (resp.^) the proper (resp. comoving) distance, a the scale 
factor and 6(r) the density contrast. At this level we clearly see 
the projection nature of weak lensing signal. This means that a 
given /e-peak can be the result of a single (5-peak (e.g. a cluster 
of galaxies) or the sum of two or more less pronounced 5-peaks 
(e.g. groups or filaments aligned along the line of sight...). 



3. The data 

3.1. Imaging and Photometry 

The CFHTLS Deep survey is intimately linked to the 
Supernovae Legacy Survey (SNLS, Astier et al. 2006) as they 
share the same data. In practice, for each observed field data are 
taken sequentially every 3-4 nights, 6 months per year and in 
4 observing bands (g', r' , i',z'). Additional u* data are included 
but they are not part of the SNLS and do not require time sam- 
pling. Most of the data have a seeing requirement limited to 
0.9". Data acquisition is still under progress at CFHT so the 
depth of the Deep fields is still improving. Data processing (as- 
trometry, photometric calibration, final stacking and production 
of catalogues) is performed at Terapix 2 for the CFHTLS com- 
munity. The final data are released regularly by the CADC 3 . 
The present analysis is based on 2 sets of data released subse- 
quently (the details of the release contents can be found on the 
Terapix Web site). For the weak lensing signal extraction we 
used the T0002 data while the deeper images and catalogues 
from the T0003 release were included for the photometric red- 
shift production. The shear analysis (see below) was already 
performed on the T0002 data release at the time of T0003 data 
release. Since the band images did not gain much depth be- 
tween T0002 and T0003, the weak lensing analysis, which is 
based on shape measurements of galaxies and is not sensi- 
tive to magnitude completeness, remained almost unchanged. 
Conversely the gain in exposure time and limiting magnitudes 
for other filters was extremely useful for photometric redshifts 
(see Sect. 3.3). 

The Deep survey is made of four independent patches 
called Dl, D2, D3 and D4. For each field and filter, 
Table 1 summarises the main observational properties of the 
T0002/T0003 release data in terms of coordinates, seeing, ex- 
posure time and depth. Because of the presence of bright stars, 
fields boundaries, defects in the CCDs and gaps between them 
with low signal-to-noise ratios, a substantial part of the images 
cannot be used for weak lensing analysis. Hence, the masked 
regions generally result in a loss of 20% of the field area. In 
Fig. 1, holes in the distribution of stars reveal the masked re- 
gions. The effective usable area is given in Table 1. The total 
working area for the weak lensing analysis is 3.61 deg 2 . 

3.2. Shear measurement 

The coherent stretching of background sources produced by the 
weak lensing effect is measured using i' band images. From 
one field to another the depth of the catalogues is slightly vary- 
ing but for better coherence between the fields, we define a 
common magnitude cut for the selection of background sources 
22 < < 26. In addition, close galaxy pairs with angular sepa- 
rations less than 3" are discarded to avoid blended systems with 
biased ellipticity measurements. Small galaxies with a half flux 
radius r/, smaller than that of stars are also rejected. 

The reliability of shear measurements is expected to be 
comparable to the current cosmic shear survey analyses like 



2 http://terapix.iap.fr/ 

3 http://cadcwww.hia.nrc.ca/ 
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Table 1. Summary of CFHTLS Deep data, release T0003 (T0002/T0003 for i' band). The limiting magnitudes correspond to 
50% completeness and are expressed in AB system. Seeing is measured in T0002 release i' images. 





Dl 


D2 


D3 


D4 


a 2000 


02 h 25 m 59 s 


10 h 00 m 28 s 


14 h 19 m 27 s 


22 h 15 m 31 s 


^2000 


-04° 29' 40" 


+02° 12' 30" 


+52° 40' 56" 


-17° 43' 56" 


zf 


25.0 


24.9 


25.1 


25.0 


i' 


25.9/25.9 


25.4/25.7 


25.9/26.2 


25.7/26.0 


r 1 


25.0 


26.0 


26.4 


26.4 


g 1 


26.4 


26.2 


26.6 


26.3 


u* 


26.5 


26.1 


25.9 


26.5 


i' seeing (") 


0.91 


0.95 


0.91 


0.87 


Area (deg 2 ) 


0.93 


0.89 


0.92 


0.86 



Dl D2 D3 D4 




Fig. 1. Stellar ellipticities in the CFHTLS Deep fields for Dl... D4 from left to right. Upper row: polarisation field before and after PSF 
anisotropy correction (respectively red upper and lower black panels.). Lower row: projection of the stellar ellipticities in the (ei,e2) plane 
before and after PSF anisotropy correction (black crosses and red dots respectively). 



the one already performed by Semboloni et al. (2006) and 
Hoekstra et al. (2006) who found similar results with inde- 
pendent pipelines in the same fields. Throughout this paper, 
we report results of shear analysis performed in the i' band 
which presents the best balance between depth and seeing. 
However we checked that we can extract a similar signal in the 
other noisier g',r' and z' bands. This has also been assessed 
by Semboloni et al. (2006) who report a similar cosmic shear 
signal in the deep fields using both r' and i' bands. 



Blurring and distortion of stars and galaxies produced by 
instrument defects, optical aberration, telescope guiding, atmo- 
spheric seeing and differential refraction are corrected using 
the Point Spread Function (PSF) of stars over the whole field. 
Several correction techniques and control of systematic errors 
have been proposed over the past 10 years (see e.g. Mellier 
1999; Bartelmann & Schneider 2001). In the following we use 
the most popular KSB method initially proposed by Kaiser 
et al. (1995). Several teams have already demonstrated that this 
technique can correct systematics down to a lower value than 
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the very weak cosmic shear signal (van Waerbeke et al. 2005; 
Heymans et al. 2006). It is therefore well suited for this analy- 
sis too. The KSB implementation used here is identical to that 
of Gavazzi et al. (2004). 

The observed ellipticity components e ob ^ , are made of the 
intrinsic ellipticity components e^°, and linear distortion terms 
that express the instrument and atmospheric contaminations 
and the contribution of gravitational shear to the galaxy ellip- 
ticity. Each ellipticity component is transformed as: 



with PL = p± - P sm 



•T I psm 



(9a) 
(9b) 



where g is the reduced gravitational shear, P sm is the smear po- 
larisability, P sh the shear polarisability and P 7 the isotropic cir- 
cularisation contribution to the final smearing. Sources are de- 
tected with SExtractor (Bertin & Arnouts 1996) but shape pa- 
rameters are calculated with Imcat 4 . Because the noise present 
in these measured quantities is important, all these tensors are 
simplified to half their trace (see e.g. Erben et al. 2001). 

(jm) and q" are measured from field stars. Their spatial 
variation across the field is fitted by a second order polynomial, 
applied individually to each one of the 36 CCDs composing 
the MegaCam focal plane. Stars are selected in the magnitude- 
Th plane, as usual, q* is the anisotropic part of the PSF, which 
is subtracted from observed ellipticities. The residuals for stars 
are shown in Fig. 1 . After correction, these latter are consistent 
with a cr y =i 0.004 rms featureless white noise. In Sect. 4.1, 
we present mass reconstructions inferred from the shear field 
measured on distant source galaxies. If we perform the same 
reconstructions on stars, we only get white noise as expected 
from a correct PSF smearing correction. Its rms is cr K ~ 0.001 
when smoothed with a Gaussian kernel of width 1 arcmin, 
which is much below the signal we are interested in (typically 
0.01 < k < 0.5, see below). 

The smearing part of the PSF contained in the P 7 term de- 
pends on the magnitude of the object and on its size as com- 
pared to the seeing disk. To optimally extract P y , we derived 
it from an averaged value over its 40 nearest neighbours in 
the magnitude-r/, plane. The variance of ellipticities inside this 
neighbourhood is then used as a weighting scheme for the shear 
analysis. The weight assigned to each galaxy is : 



1 



py2 



W: = 



(10) 



where ctq prevents from over-weighting some objects. It is set 
to the ID intrinsic dispersion in galaxy ellipticities cr e = 0.23 
and cr, is the observed dispersion of ellipticities over the 40 
nearest neighbours in the magnitude-r/, plane. 

At this level, we have 132,000 (resp. 104,000, 162,000 and 
114,000) galaxies with reliable shape parameters in the Dl 
(resp. D2, D3 and D4) field leading to a source surface num- 
ber density of rt\, g = 38.0 (resp. 30.6, 35.4, 34.5) arcmin -2 . 
These values are much higher than the usual ones in weak lens- 
ing studies, which turn around 15-20 galaxies arcmin -2 . The 



magnitude cut i' < 26 explains these high densities although 
it makes difficult an accurate determination of the redshift dis- 
tribution of such faint objects. In the following we shall re- 
fer to this source catalogue as CA. Because several galaxies 
have large uncertainties on their shape parameters, one should 
correct this density by considering the effective density which 
would have been achieved if all uncertainties were limited by 
the intrinsic dispersion in source ellipticity. More precisely if 
we define A^ff = 2 iiCe/Cej) 2 , the effective source number den- 
sity would then be «£ ff = 25.3 (resp. 20.7, 21.6, 21.0) arcmin -2 . 
This quantity represents what could be achieved under deep 
space-based observing conditions. 

Fig. 2 illustrates the effect of our weighting scheme (Eq. 
10) in the magnitude distribution of sources. The weighting 
scheme efficiently down-weights objects fainter than i' ~ 24 
and most of the signal is carried by galaxies brighter than 
i' ~ 25. Actually i' = 24 is the magnitude above which shape 
measurement errors start to increase above the intrinsic scatter 
of source ellipticities. 



10 5 



10 4 



io- 



10' 




18 20 22 24 

i' band magnitude 



28 



Fig. 2. V band magnitude distribution of background sources. For 
clarity we represent the whole magnitude distribution although the 
CA catalogue only contains galaxies in the range 22 < i' < 26 . The 
dashed red line illustrates the effective change in magnitude distribu- 
tion of lensed sources when accounting for the weighting scheme of 
Eq. (10). The weights are normalised so as to conserve the total num- 
ber of objects. 



3.3. Photometric redshifts 

The number of filters available in the CFHTLS Deep images al- 
lows a direct estimate of the redshift of each source according 
to its Spectral Energy Distribution (SED) from u'g'r'i'z' bands. 
We used the publicly available catalogue of photometric red- 
shifts recently carried out by Ilbert et al. (2006) for which the 
accuracy has been extensively assessed against spectroscopic 
data 5 . In some cases Ilbert et al. used additional non-CFHTLS 
data in other photometric bands to improve the precision of 



http : //www. ifa.hawaii . edu/~kaiser/imcat/ 



5 http://terapix.iap.fr/article.php?id_article=586 or 
http : / /cencos . oamp . fr/cencos/CFHTLS/ 
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their photometric redshifts. We refer the reader to this work 
for further accounts. Here we briefly mention that photometric 
redshifts have been checked to be reliable enough for our pur- 
pose down to i' ~ 25. The authors mention that the accuracy is 
reduced at redshift Zphot > 1.5. 

In the following we use Ilbert et al. (2006) photometric red- 
shifts in three distinct ways, each of them requiring different 
precision on z p h t- 

- To estimate the redshift distribution of the whole sample 
of background galaxies CA, we use photometric redshifts 
down to i' < 26 which is a rather faint limit for accurate 
photometric redshifts. However the magnitude-dependent 
weigthing scheme of Eq. (10) means that the net contribu- 
tion of lensed sources peaks at if ~ 24 and drops quickly for 
fainter objects (see Fig. 2). Fig. 3 shows that the weighted 
redshift distribution of i' < 26 sources is close to that of the 
unweighted distribution of i' < 24 galaxies for which pho- 
tometric redshifts are accurate. In addition we see on this 
figure that 80% sources are at redshift below 1 .5. Therefore 
the relatively degraded precision of photometric redshifts 
beyond z = 1.5 is not an important concern for the de- 
termination of the redshift distribution of CA. The mean 
weighted redshift of sources is = 0.92 (whereas it is 

- 1.01 if weights are neglected). 

- In Sect. 5.1 we will need to know the redshift of back- 
ground sources individually. We discuss there the effect of 
changing the limiting magnitude or redshift of our sample. 

- In Sect. 5.2, we identify the optical counterpart of structures 
detected by weak lensing. We use the photometric redshift 
distribution of the bright galaxies present in the structure as 
a proxy for the structure redshift itself. In this case we will 
consider bright if < 23 lens galaxies for which photometric 
redshift are very accurate <x. /( 1 + z) ~ 0.03 down to z ~ 1.5 
(Ilbert et al. 2006). 



4. Mass reconstructions 

4.1. Convergence maps from observed shear 

From the source catalogue CA defined in Sect. 3, we can infer 
the shear field y(6) and deduce the associated convergence field 
k(6). They are related by: 

k(0)= f K(0-&)*y(&)d 2 $, (11) 

Jr 2 



where K(ff) 



n (6 



2 is a complex convolution kernel (Kaiser 



& Squires 1993) (hereafter KS93). The shear field is smoothed 
with a Gaussian filter G{6) oc exp(-^r) with 6 S = 1 arcmin. 
The convergence field is consequently smoothed by the same 
filter. The resulting convergence maps present correlated noise 
properties (van Waerbeke 2000). 

Annuel \ 40 2 , 



(K n {$)K n (& + 0)) 



cr^Ann^fr]) characterises the noise level. We measured a 
value 0.0196, 0.0225, 0.0202 and 0.0221 for Dl, D2, D3 and 
D4, respectively. 




3.5 



Fig. 3. Redshift distribution of sources brighter than i' = 26. Thin 
Solid black: (resp. thick solid red) the weighting scheme in Eq. (10) 
which reduces the contribution of faint-distant sources is (resp. not) 
taken into account. For comparison the unweighted redshift distribu- 
tion of i' < 24 sources is shown (dashed blue). This latter distribution 
is renormalised to the same total number of objects as the i' < 26 
sample. 



In principle, the convergence computed from Eq. (11) must 
be real and its imaginary component should only be due to 
noise and possible residual systematics. We checked this as- 
sumption by rotating the shear field by 45° and found the re- 
constructed maps to be consistent with noise as described by 
Eq. (12). 

The KS93 inversion in Eq. (1 1) is done by a direct summa- 
tion over all sources without pixelling, smoothing and Fourier 
transforming the data. This reduces boundary and mask effects 
on mass reconstructions. Several techniques have been pro- 
posed so far since the original KS93 method. Most of them are 
useful in high shear regions (where g < 1) and for small fields 
of view. However the wide MegaCam images and the complex 
field geometry imposed by the masks make difficult, time con- 
suming and unnecessary the implementation of more complex 
techniques. In addition, van Waerbeke (2000) has shown that 
noise properties of KS93 method are well controlled and con- 
sistent with Eq. (12). 

Fig. 4 shows the convergence maps for Dl, D2, D3 and 
D4 deduced from the catalogue CA. Contours in units of the 
signal-to-noise ratio (SNR or v) are overlaid, with v defined as 



SNR 



«bg 



el 



(13) 



In the present data we detect ~ 46 positive peaks with v > 3 and 
5 peaks with v > 4. In order to avoid too much contamination 
by noise peaks but to detect as much true peaks as possible, we 
therefore fix the threshold at v = 3.5. The 14 peaks detected 
within this limit will constitute our working sample in the rest 
of the paper. We discuss in more detail the statistics of these 
peaks in Sect. 4.2 and their possible association to galaxy clus- 
ters in Sect. 5. 
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Fig. 4. Convergence maps inferred from the ellipticity field of background sources for Dl (top left), D2 (top right), D3 (bottom left) and D4 
(bottom right). The shear field is estimated by smoothing the ellipticity field of sources galaxies selected in catalogue CA. Then, Eq.(l 1) is used 
the convert y into the k field. The largest masked regions are visible as fuzzy k ~ regions. Contours levels start at 3<x with a 0.5<r arithmetic 
increase. The Gaussian filtering scale is 1 arcmin. The 14 peaks with v > 3.5 are labeled. 



4.2. Statistics of peaks 

Several authors investigated the possibility to use convergence 
peaks as clusters of galaxies candidates. Simplified analytical 
calculations based on the Halo Mass Function (as inferred from 
the Press-Schechter formalism for instance) provided the first 
predictions for wide field imaging surveys (Schneider 1996; 
Kruse & Schneider 1999). Then, thanks to the development of 
numerical simulations, quantitative estimates of projection ef- 
fects and cluster selection functions (in terms of mass and red- 



shift) became available (Reblinsky & Bartelmann 1999; Jain 
& van Waerbeke 2000; White et al. 2002; Padmanabhan et al. 
2003; Hamana et al. 2004; Hennawi & Spergel 2005; Tang & 
Fan 2005). 

The practical implementation of a Weak Lensing Cluster 
Survey (WLCS) requires the control of noise present in obser- 
vations, either due to the intrinsic ellipticity of sources or to the 
intervening large scale structures (LSS) along the line of sight. 
Although the "compensated aperture mass filter" has early been 
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proposed as an efficient filter for peak statistics (Schneider 
1996; Kruse & Schneider 1999; Schirmer et al. 2003, 2004; 
Hetterscheidt et al. 2005; Schirmer et al. 2006), it has been 
shown that such a filter may not be as efficient as an optimised 
filter which would account for the contribution of LSS to the 
noise budget and the shape of the dark matter halos we aim 
at detecting (Hennawi & Spergel 2005; Maturi et al. 2005). It 
turns out that a simple Gaussian filter of width S ~ 1-2 arcmin 
is close to the optimal linear filter and has been extensively 
studied in simulations (White et al. 2002; Hamana et al. 2004; 
Tang & Fan 2005). In addition a promising multiscale wavelets 
technique has also been proposed recently (Starck et al. 2006) 
but has not been applied to real data yet. Therefore we shall 
use a simple Gaussian filter with scale 8 S = 1 arcmin as already 
applied onto the mass reconstructions of Sect. 4. 1 . 

Fig. 5 shows the cumulative number of maxima peaks 
N(> v) as well as the symmetric number of minima peaks for 
the four Deep fields. The latter curve is flipped (v — » -v) for 
an easier comparison. The net excess of maxima with v > 4 as 
compared to minima at the corresponding negative threshold 
is visible, thus showing the non-Gaussian nature of the con- 
vergence field (see also Miyazaki et al. 2002). The statistical 
significance of this excess is still low due to the large cosmic 
variance in such a small sky coverage. In addition it should be 
kept in mind that CFHTLS Deep fields of view were chosen to 
be free of any known massive nearby cluster. 

The theoretical analysis of Hamana et al. (2004) is well 
suited for a direct comparison with our results since the sur- 
vey area, the smoothing scale, and the noise properties are the 
same. We found a satisfying agreement when considering their 
Fig. 7 although the sample variance is large. 
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Fig. 5. Cumulative counts of N(> v) maxima peaks (solid black) 
and corresponding counts of N(< -v) minima peaks (dashed blue) 
curve per square degree in the Deep survey. The surrounding dotted 
lines represent lcr Poisson errors. The excess of positive maxima is 
marginally seen as compared to negative minima, thus showing the 
non-Gaussianity of the convergence field. 

An extensive study of /c-peaks statistics in the Wide 
CFHTLS survey would provide valuable cosmological infor- 
mation for cosmic shear studies and would help breaking some 



degeneracies (mainly between Q. m and cr 8 ) present in the shear 
2-point correlation function. Such an analysis is beyond the 
scope of the present work but the full implementation of the k- 
peaks statistics in the presently released CFHTLS-Wide data is 
in progress. It is noteworthy that in order to extract the associ- 
ated cosmological signal, there is no need to measure the mass 
or the redshift of each peak, neither to identify them with clus- 
ters (or projected groups, etc.). However, instead of a blind 
counting exercise of convergence peaks, we propose in the fol- 
lowing to characterise the structures responsible for the highest 
convergence maxima peaks. 

5. Properties of v > 3.5 peaks 

In this section we attempt to estimate the redshift, mass and 
luminosity of structure(s) responsible of the 14 /e-peaks with 
v > 3.5. This significance threshold is rather low so we expect 
a substantial amount of contamination by noise fluctuations. 
We focus on the physical properties of the detected peaks and 
test whether they can be associated to galaxies over-densities. 
The first step in this peak identification is to estimate a redshift, 
taking advantage of photometric redshifts. Two methods are ex- 
plored, one which is directly related to the shear signal depen- 
dency on lens redshift and which does not require an explicit 
identification of the peak with galaxies, and the other which is 
related to the identification of a localised over-density in the 
photometric redshift distribution of galaxies. Results from the 
two approaches are summarised in Table 2. 

5.1. Lens tomography 

The basics of lens tomography are the following: in the case of 
a real deflector at redshift zu the shear signal must increase in a 
characteristic way as a function of the source redshift according 
to the w(z\, z s ) = D\ s /D os term. Early applications can be found 
in Wittman et al. (2001, 2003) and Hennawi & Spergel (2005). 
The shape of the shear increase versus the source redshift al- 
lows to derive the lens redshift. In the following, we apply this 
technique to estimate the lens redshift of each peak and check 
whether the shear behaviour around peaks is consistent with 
the lensing hypothesis or is rather a noise fluctuation. 

We measure the shear profile around each peak between 
1 and 5 arc minutes from the centre, using sources having 
i' < 26 and their individual photometric redshifts. Assuming 
that the lens mass distribution at redshift z\ follows a Singular 
Isothermal Sphere (SIS) profile, the shear is simply 



y(6»,z s ) = w(zi,z s ) 



29 ' 



(14) 



The Einstein radius 9e is related to the characteristic cluster 
velocity dispersion by <r v = 186.2kms~ (0e/1")^ 2 - We fit the 
shear function for the unknown lens redshift zi and Einstein 
radius 6 E by minimising ax 2 (zi, 6 E ) of the form 



(15) 



where w; = w(zi,Z s ,d, o~e,i is given by Eq. (10) and e,j is the 
tangential component of ellipticity. The dependency of x 2 with 
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Fig. 6. Increase of shear signal strength (i.e. equivalent Einstein radius) as a function of source redshift around peak Cl-02. The source catalogue 
is split into 10 quantiles. Contours show 1, 2 and 3 <x CL regions around the best fit 8e and zi. Peak Cl-02 exhibits the expected profile for a 
cluster with velocity dispersion <r v = 760 ± 1 10 km s -1 and redshift z\ = 0.52+jJ ]| with^ 2 /dof = 0.6 (thick line). The left panel uses photometric 
redshifts of sources down to i' = 26 whereas the right panel is limited by i' < 24.5, showing that including faint sources does not affect our 
analysis in a statistically significant way. 



6>e can be easily removed by considering that ^ = for any zi 
if 9e satisfies 

— ^6^. 

9e = ^ (16) 

y. J!L- 

It can be inserted into Eq. (15) to give 

2, , V j|] m 

x (zi) = h~2 — • (I?) 



Y-4 



a 2 ., 



In order to illustrate the method we plot on the left panel 
of Fig. 6 the value of the Einstein radius measured in the 10 
z s quantiles of sources between 1 and 5 arcmin from the cen- 
tre of the peak Cl-02, detected with v = 5.5. The increase of 
6>e with redshift z s is clear and allows an unambiguous iden- 
tification of the lens redshift. Error bars come from the scatter 
in observed (measurement+intrinsic) ellipticities as determined 
by Eq. (10). Contours show the 68.3%, 95.4% and 99.3% CL 
regions for z\ and 9e- The thick curve represents the function 
&e x w(z\,Zs) for the best fit values z\ = and 6 E = 

17^5 »" leading to a velocity dispersion cr v = 760 + 1 10 km s . 
Because of the rather large statistical uncertainties, the detailed 
radial shape of the shear profile introduced in the x 1 fit is 
not much important. For significant peaks in the reconstructed 
mass maps a redshift estimate directly evaluated from the lens- 
ing properties of the peaks is a viable method. It could in par- 
ticular be applied to any putative "dark clump" where no bright 
galaxies can be associated with the mass peaks. 

The right panel of Fig. 6 shows that limiting the analysis 
to brighter sources (including i' < 24.5 photometric redshifts) 
does not make significant changes. We checked that this is also 
true for other peaks. This suggests that, even if photometric 
redshift of i' > 24.5 galaxies may suffer larger uncertainties, 
the effect on lens tomography, and lensing in general, is negli- 
gible. This is due to the saturation of the curve Di s /D os at high 



z s and to the weighting scheme (Eq. 10), as already mentioned 
in Sect. 3.3. Table 2 presents the constraints on cr v and z\ given 
by tomography for the 14 peaks. Most of them have fitted ve- 
locity dispersion values limited to cr v < 700 km s . There is no 
massive cluster with cr v > 800kms~' below redshift ~ 0.7 in 
the Deep survey. This is not a surprise as the Deep fields were 
initially selected for their lack of well identified Abell clusters 
for example. In addition we observed that if the signal-to-noise 
ratio is too low or tomography provides a bad,y 2 /dof > 2, then 
the inferred lens redshift is systematically z\ — 0. Two peaks 
(Cl-04 and CI- 12) are such that tomography fails in constrain- 
ing lens redshift and velocity dispersion. However this is may 
also be attributed to the fact that these two peaks are noise fluc- 
tuations rather than being produced by real galaxy clusters. 

5.2. Optical counterparts 

Here again we use photometric redshifts to check whether an 
optical counterpart can be assigned to our 14 high convergence 
peaks with v > 3.5. 

We first examine galaxies in a circular aperture of 2 arcmin 
around each peak. This radius corresponds to a linear physical 
scale of 400 to 860 kpc for a lens redshift ranging from 0.2 
to 0.7 respectively. Therefore this is the radius within which 
the highest density of bright galaxies is expected, most of them 
being early-type galaxies. In addition, most of these galaxies 
are localised in a narrow range of photometric redshifts once 
the background z p h t -distribution is subtracted. The background 
is defined in the region beyond 6 arcmin of all peaks. In this 
preliminary step we only consider galaxies brighter than i' = 
23 and with a best fit SED template type T of type earlier than 
the spectral type T = 44 6 . 



6 The best-fit galaxy model T given by Ilbert et al. (2006) there- 
fore excludes starbusts and Irregular galaxies for which photometric 
redshifts are less reliable. 
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Fig. 7. Photometric redshift distribution of galaxies around convergence peaks with a candidate optical counterpart. The green dashed curve is 
the lens redshift as constrained by tomography (Sect. 5.1), showing a remarquable agreement with z phot excesses for the most significant peaks. 
The total number of galaxies in the peaks should not be seen as richness since we did not correct for the presence of masks and the increasing 
incompleteness as z increases. 



9 peaks out of 14 meet this criterion whereas there is no 
clear optical associable counterpart for the remaining 5 peaks 
(rows filled with dashes in Table 2). For each of the former 9 
/c-peaks, we define the cluster redshift as the location of the 
most prominent z p h t excess peak. We iteratively apply a 4<x 
clipping to remove outliers and get a reliable redshift and its 
corresponding statistical uncertainty. Distributions can be seen 
in Fig. 7. The tomographic z\ probability distribution function 
of Sect. 5.1 (marginalised over velocity dispersion) is overlayed 
for comparison. We see a remarquable agreement for the most 
significant peaks. 

The peculiar case of CI- 14 deserves a special attention. Cl- 
14 clearly exhibits a bimodal distribution that lens tomography 
is unable to reveal, probably because of the low detection level. 
The first excess is at z — 0.15, and is d — 75" away from 
the convergence peak location. It is consistent with an XMM- 
LSS detection at a similar redshift (see Sect. 5.4). However the 
second peak, at redshift z — 0.57, is much closer to the conver- 
gence peaks (d = 14") although it does not match an XMM- 
LSS peak. It is likely that both components may contribute to 
the overall convergence. This is a clear example of the pro- 



jection effects already mentioned at the end of Sect. 2. In the 
following we will consider separately these two components, 
labeled Cl-14a for the z 0.15 peak and Cl-14b for the fur- 
thest one. 

We also define the luminosity- weighted optical centre using 
the bright galaxies (/' < 23) in the z p h t excess range. Optical 
centres are less than one arcmin away from the convergence 
peak location, as expected from the spatial resolution of the 
mass maps. The observed offsets are reported in column A(a, 5) 
of Table 2. From the new defined centre and the present cluster 
photometric redshift, we re-estimate the lensing velocity dis- 
persion (central part of Table 2). The velocity dispersions of Cl- 
14a and b are fitted simultaneously. Although they are slightly 
correlated, we can deblend the system and put constraints on 
both components. In any case, CI- 14b is the dominant contri- 
bution to the convergence peak. 

In order to relate lensing velocity dispersion and the excess 
of foreground galaxies we measure the number A of M r < Mo 
galaxies in a fiducial physical radius Ro = 1 .43 Mpc and a slice 
fephot = ±0.10, with Mo - -20 in the rest-frame r' band. 
The most distant of our peaks being at z ~ 0.74, the sample 
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Notes: (a) ROSAT cluster (22 h 16 m 56.2 s , -17° 25' 25.5") at z = 0.13 (de Grandi et al. 1999), (b) XMM-LSS cluster at z = 0.31 (XLSSC13), 
(c) XMM-LSS cluster at z = 0.26 (XLSSC25), (d) XMM-LSS cluster at z = 0.29 (XLSSC22), (e) XMM-LSS cluster at z = 0.14 (XLSSC41). 
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optical counterpart. Rows filled with "X" are in the Dl field part of the XMM-LSS survey but not detected in X-rays. Cl-02 is part of Dl but 
lies in a region lost by the XMM-LSS survey (due to high flare rates). The second component of Cl-14 is not detected by the XMM-LSS as it 
may be hidden by the foreground z = 0.14 component. 



of photometric redshifts that will be close to complete at Mq 
corresponds to i' ~ 23.5, which is a conservative limit for ac- 
curate photometric redshifts. Their contribution to the L > L* 
luminosity function is very low anyway. Assuming a luminos- 
ity function with slope a, the total luminosity in the aperture 
is 

Laper = LA ^ + Qr, °^ 1.24AL., (18) 
L( 1 + a, q) 

with q = lO 04 ^--^) and assuming a = -1.09, AT, = -21.21 
(Blanton et al. 2003). This procedure limits the impact of bad 
photometric redshifts as compared to a simple summation of 
galaxies luminosities. Because photometric redshifts are im- 
perfect, we need to subtract a background contribution which is 
evaluated using galaxies outside 1 .5Rq of all peaks in the cor- 
responding redshift slice. We also compensate for the aperture 
area lost in masked regions. CI- 14a and CI- 14b are considered 
independently. 

The correlation between lensing velocity dispersion and 
richness is shown in Fig. 8 for the eight systems Cl-01, Cl-02, 
Cl-03, Cl-05, Cl-06, Cl-07, Cl-08, Cl-12 and the Cl-14 which 
is split into its two componants. We also add a peak with snr 
v = 3.4 as it is part of the XMM-LSS cluster sample (see Sect. 
5.4). The range of richness and velocity dispersion probed by 
our sample is quite narrow so it is difficult to put tight con- 
straints on the scaling relation between these quantities. A lin- 
ear regression yields : cr v ~ 251 km s _1 A - 21 . The slope is how- 
ever poorly constrained. 

For our assumed isothermal mass profile, the projected 
mass enclosed in radius Rq is: 

M(< R Q ) = nc^Ro/G . (19) 



From the data we find <r = 96 + 7 km s~' A 1 ^ 2 (solid green line 
of Fig. 8). Combining Eq. (18) and Eq. (19), this translates into 
a r' band mass-to-light ratio M/L r = 168+24 /! 7 o(M/L) G inside 
the aperture Rq. The z = 0.74 peak Cl-08 slightly departs from 
this relation but this could partly be due evolution of galax- 
ies luminosities. In addition, weak lensing mass estimates for 
such high redshift clusters are particularly sensitive to the high 
redshift sources. Since the accuracy of photometric redshifts is 
more hazardous at z s > 1 .5 the lensing-inferred velocity disper- 
sion of this cluster is more uncertain. Otherwise our M/L ratio 
is consistent with recent estimates (see e.g. Yee & Ellingson 
2003; Popesso et al. 2006; Bardeau et al. 2005). 

5.3. Merits of lens tomography 

It is instructive to compare the redshift arising from lens to- 
mography (ztomo) and the redshift z p hot given by an excess of 
foreground galaxies (i.e. cluster members identified by their 
photometric redshift). It is not a surprise that most WLCSs de- 
tections peak at <z to mo) - (z P hot) - 0.25 because is corresponds 
to the redshift range in which gravitational lensing is most effi- 
cient for the source population we are considering. 

In Fig. 9 we plot the comparison between these two red- 
shift estimates for the 8 clusters having a well identified optical 
counterpart (namely clusters Cl-01, Cl-02, Cl-03, Cl-05, Cl- 
06 Cl-07, Cl-08, Cl-12). Cl-14 was discarded due to its appar- 
ent complexity. Although the statistic is quite small and errors 
on Ztomo are large, there is a remarquable agreement between 
both redshifts for the high signal-to-noise systems, which is 
very encouraging. In addition we already pointed out that Cl- 
04, for which tomography fails in giving a lens redshift, has no 
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Fig. 8. Correlation between richness and lensing velocity dispersion 
for Cl-01, Cl-03, Cl-05, Cl-06, Cl-07, Cl-08, Cl-12, Cl-14a-b and 
XLSSC44. The dotted red line is a raw linear fit leading to the scaling 
o"„ ~ 251 km s~ 1 A 2 1 . If we assume a constant mass-to-light ratio, we 
find <r v = (95 ± 7)A 1/2 km s _1 i.e. M/L r = 168 + 24 h 10 (M/L) G (solid 
green line). 
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Fig. 9. Comparison between lens tomography and direct photometric 
redshift methods to estimate the redshift of the v > 3.5 /(-peaks. Only 
those peaks with a clear optical counterpart are included in the figure. 
Red diamonds represent the 4 most significant v > 4 peaks whereas 
Blue squares code for the remaining 4 lower snr peaks. Due to its 
complexity CI- 14 is not included in this figure. 



optical (nor X-ray) counterpart and is likely a false detection. 
Tomography helps eliminating such cases. 

In order to improve this correlation, it will be important to 
improve the quality of photometric redshifts, especially for the 
faintest background sources. It is not clear how much the catas- 
trophic redshifts in the source catalogue perturb lens tomogra- 
phy, but certainly for accurate tomographic redshift estimates 
it will be important to increase the number of filters for photo- 
metric redshifts, especially in infrared bands (Bolzonella et al. 
2000). In addition, it is noteworthy that our sample is made of 
low mass clusters (cr v ~ 400 - 600kms~'). For more massive 
clusters, like those expected in the CFHTLS-Wide survey, the 
method will greatly gain in accuracy and reliability (Hennawi 
& Spergel 2005). 

5.4. Comparison to X-ray data 

The CFHTLS -Deep Dl field is part of the XMM-LSS survey 7 . 
We therefore took advantage of the publicly available X-ray 
database to cross-correlate our sample of K-peaks/clusters with 
those X-ray detections published in (Valtchanov et al. 2004; 
Willis et al. 2005; Pierre et al. 2006). 

The matching is very good. Over the seven v > 3.5 peaks 
detected in Dl, four are XMM-LSS clusters with luminosity 
1.5 10 43 < L x ,boi < 6.5 10 43 erg/s and temperature 1 < T x < 2 
keV. However the most pronounced v = 5.5 Dl peak, namely 
Cl-02, is not part of the XMM-LSS sample because it falls in 
a region lost by the X-ray survey (pointing G12, Pierre et al. 
2006). Peaks Cl-04 and Cl-10 do not exhibit any optical coun- 
terpart and are not detected in X-ray either. This confirms that 



7 http : //vela . astro . ulg . ac .be/themes/spatial/ 

xmm/LSS/index_e .html, see also http://13sdb.in2p3.fr: 
8080/13sdb/ 



Cl-04 and Cl-10 are likely false detections due to noise fluctu- 
ations. 

There are 9 publicly available XMM-LSS clusters in the Dl 
field of view in the classes C1-C2 8 . 4 of them are part of our 
weak lensing cluster sample although we note that a v = 3.4 
peak at a = 02 h 24 m 31.7 s and 5 = -04° 13' 55" is also part of 
the XMM-LSS sample (XLSSC44 at z = 0.26, T x = 1.37^ 
keV). This cluster has been missed by the weak lensing sur- 
vey because it does not meet the v > 3.5 detection threshold. 
Its lensing velocity dispersion is <T v j e „ s = 466^[^kms _1 . The 
remaining 4 XMM-LSS clusters which are not part of our sam- 
ple: XLSSC38 (z = 0.58, T x unknown), XLSSC11 (z = 0.05, 
T x = 0.64+° ii keV), XLSSC29 (z = 1.05, T x = 4.07^-™ keV) 
and XLSSC5 (z = 1.05, T x = 3.67+ 3 |0 keV ) are either very 
low or very high redshift clusters for which lensing efficiency 
is low. Therefore it is not surprising that they are missing in our 
weak lensing sample. 

Note also that the peak Cl-01, located in D4 and which is 
the strongest peak matches an X-ray detected ROS AT cluster at 
redshift z = 0.13 and luminosity Lx,[0.5-2.o] = 14.5 + 2.5 x 10 43 
erg s _1 . (de Grandi et al. 1999). 

We now compare lensing velocity dispersion and X-ray 
temperature for the five Dl XMM-LSS clusters Cl-03, Cl- 
05, Cl-07, CI- 14 and XLSSC44. More precisely and following 
Cypriano et al. (2004), we compare the lensing velocity disper- 
sion inferred temperature ksT\ ens = fjm.H<T$ of dark matter par- 
ticles to the X-ray temperature of hot ICM gas particles, ymu is 
the mean particule weight per free electron. It is set to the same 
value yu = 0.61 as Cypriano et al. (2004). We also include in 
the comparison data from the study of Bardeau et al. (2006) of 
1 1 X-ray luminous clusters at z = 0.2. Under the assumption of 
energy equipartition these temperatures would be equal. If non 

8 4 additional C3-class clusters could be added to this list but their 
detection and physical properties are more hazardous 
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Fig. 10. Correlation between X-ray and lensing-inlerred temperature 
k B T icm = fjm H crl (black diamonds) for Cl-03, Cl-05, Cl-07, CI- 14 and 
XLSSC44. The results of Cypriano et al. (2004) on massive clusters 
are reported for comparison (red stars) as well as those of Bardeau 
et al. (2006) (green crosses). The bisectrix line (dotted) represents en- 
ergy equipartition between hot ICM gas and dark matter + galaxies 
collisionless fluid. 

gravitational sources of gas heating/cooling are at work we ex- 
pect some departure from this relation Tx + Ti ens . Conversely 
the mass (and thus cr v and T\ eni ) of shear-selected clusters may 
be increased by projections of unrelated material along the line 
of sight. 

Although the statistics is rather poor, Fig. 10 suggests that 
shear-selected clusters are well aligned onto the bisectrix Tx = 
7i ens . Gas and collisionless particles share the same amount of 
energy. This behaviour seems to be less true for massive clus- 
ters. For example Cypriano et al. (2004) found that for Tx > 8 
keV, the gas is hotter than expected by pure gravitational ef- 
fects. This supports the presence of off-equilibrium physical 
processes (unrelaxed clusters, merging). 

The low redshift component of CI- 14 (CI- 14a) has been 
used for Fig. 10 although it should be considered with caution. 
Both lensing and X-ray properties of this cluster may be con- 
taminated by projection effects: CI- 14b seems to dominate the 
lensing signal whereas the nearby CI- 14a dominates the X-rays 
signal. 

6. Discussion and conclusion 

In this work we attempted to analyse the weak lensing signal 
in the 4 deg 2 images of the CFHTLS Deep fields. For a proper 
signal extraction we payed special attention to the removal of 
residual systematics and got a large sample of distant lensed 
sources, thanks to the exceptional depth of the images. We then 
used standard KS93 inversion technique to infer the projected 
surface mass density (i.e. the convergence field) and focused on 
maxima peaks with a signal-to-noise ratio v > 3.5. We found 
14 such peaks and discussed the possibility to use the statistics 
of maxima peaks as a test for the non-Gaussianity of the con- 
vergence field. The lack of massive clusters is not surprising 
since the Deep field of view is only 4 deg 2 and the 4 MegaCam 
fields were selected a priori as free from already known rich 



clusters. It turns out that half of the v > 3.5 convergence peaks 
(either with or without an optical counterpart) are in Dl and at 
redshift z ~ 0.28, i.e. in the redshift range with best lensing effi- 
ciency. How much of this excess is due to the 10% lower noise 
level in Dl relative to D2, D3, D4 or is pure sample variance 
(enhanced by the strong spatial clustering of galaxy clusters)? 
The latter is our favoured explanation since a 10% change in 
SNR for Dl peaks would not significantly change the ranking 
of peaks in Table 2. At z ~ 0.28, 1 deg scale corresponds to 
15 Mpc and the cosmic variance of the cluster-cluster cor- 
relation function is still high. 

Looking for excesses of galaxies around convergence peaks 
we found that ~ 35% of our v > 3.5 candidates turn out to 
be false or inconclusive detections. Our cluster candidates are 
not very massive systems but look more like small clusters / 
large groups having 400 < cr v < 600 km s -1 . Most of them 
lie at redshift ~ 0.3 although we found a clear detection at 
z = 0.74. All the Dl XMM-LSS clusters that lie in the lens- 
ing relevant redshift range 0.1 < z < 0.6 are detected with a 
SNR v > 3.4. The completeness of WLCSs is however lower 
than X-ray techniques for clusters detections. If one aims at re- 
ducing the amount of false detections (higher efficiency), the 
sample completeness of WLCSs turns out to be very low. This 
has already been pointed out in simulations (e.g. White et al. 
2002; Hamana et al. 2004; Hennawi & Spergel 2005). 

Projections effects were observed in at least one of our de- 
tections (CI- 14). The other systems turn out to have lensing- 
inferred mass (or velocity dispersion) properties consistent 
with their optical (M/L r ~ 170/i 7( )(M/L) o ) and X-rays (cr 2 oc 
Tims - Tx) counterparts. 

We used lens tomography around /e-peaks to estimate the 
deflector redshift as well as its velocity dispersion. To this 
end we made an extensive use of photometric redshifts in 
the sample of background sources. With accurate photometric 
redshifts, tomography can improve weak lensing cluster sur- 
veys capabilities since it helps distinguishing real clusters and 
noise fluctuations. The agreement between tomographic red- 
shifts and photometric redshifts of cluster members is remar- 
quable for the most significant peaks. 

Although the field area of the CFHTLS Deep survey is not 
wide enough for cosmological application, we have demon- 
strated that CFHTLS image quality is well suited for WLCSs 
as it is the case for cosmic shear signal extraction (Semboloni 
et al. 2006). The full implementation of this technique to the 
CFHTLS Wide survey is on-going and will provide us with 
a few hundred shear-selected clusters. The large sky coverage 
will balance the lower density of background sources as com- 
pared to the Deep fields and will clearly favour the detection 
of higher mass systems with velocity dispersion in the range 
700- ^OOkms" 1 . Given our findings for the Deep fields, we 
can forecast that a detection threshold as high as v > 4.5 will 
be required for a robust WLCS in the shallower Wide survey. 
This will be done at the expense of finding intermediate mass 
structures with cr v ~ 400 - 700kms~ 1 . 

In addition half of this survey (the Wl region) will also be 
covered by the XMM-LSS survey. This will give the necessary 
calibration of scaling relations between mass and direct observ- 
ables for clusters of galaxies to be used as efficient cosmologi- 
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cal probes. A more detailed comparison with the performances 
of ongoing other clusters survey techniques (optical, SZ with 
Planck) will also become possible. The low completeness of 
WLCSs is balanced by their well controlled selection function 
since one needs cosmological simulations with relatively low 
resolution and essentially no gas physics. 

The opposite approach is also possible. We mentioned in 
Sect. 4.2 that it would be more interesting to use the statistics 
of /c-peaks as a test of the non Gaussianity of the convergence 
field. This contains complementary information on cosmolog- 
ical parameters relative to the cosmic shear 2-point correlation 
function. Like e.g. the skewness it helps breaking the observed 
degeneracy between Q. m and cr 8 with shear correlation func- 
tions. In this respect it is not necessary to check whether indi- 
vidual peaks are real or false detections nor to assign a redshift 
and mass with expensive follow-up for each convergence peak. 

Both WLCSs (with cluster identifications) and raw /(--peaks 
statistics are complementary applications of weak gravitational 
lensing. They both will soon provide new insightful constraints 
on the evolution of large-scale structure driven by Dark Matter 
and perhaps giving important clues on the behaviour of Dark 
Energy as a function of redshift. 
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